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Abstract 



We develop a theoretical frame for the study of classical and quantum 
gravitational waves based on the properties of a nonlinear ordinary dif- 
ferential equation for a function (j{rj) of the conformal time r], called the 
auxiliary field equation. At the classical level, (j{jf) can be expressed by 
' means of two independent solutions of the "master equation" to which the 

' perturbed Einstein equations for the gravitational waves can be reduced. 

, At the quantum level, all the significant physical quantities can be 

formulated using Bogolubov transformations and the operator quadratic 
Hamiltonian corresponding to the classical version of a damped para- 
metrically excited oscillator where the varying mass is replaced by the 



^ ' square cosmological scale factor 0^(77). A quantum approach to the gen- 



eration of gravitational waves is proposed on the grounds of the previous 
77— dependent Hamiltonian. An estimate in terms of (j(ri) and a(ri) of the 
destruction of quantum coherence due to the gravitational evolution and 
an exact expression for the phase of a gravitational wave corresponding 
' to any value of 77 are also obtained. We conclude by discussing a few 

applications to quasi-de Sitter and standard de Sitter scenarios. 
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1 Introduction 



In a notable paper by Grishchuk and Sidorov ^ , it was shown that reUc gravi- 
tons can be created from the vacuum quantum fluctuations of the gravitational 
field during cosmological expansion and can be interpreted as squeezed quantum 
states of the gravitational field, in analogy to what happens for squeezed states 
in quantum optics |2I4| . A systematic treatment of the particle creation mech- 
anism is contained in jS]. The theory of particle creation is essentially based 
on the Bogolubov transformations, exploited primary in the formulation of the 
process of squeezing for the electromagnetic field. An important aspect of the 
theory of particle creation is the possible loss of coherence in quantum grav- 
itational theories. In a sense, the exploration of decoherence in gravitational 
theories has been more concerned with the quantum (gravity) effects which can 
manifest themselves even at scales others than the Planck one. The growing 
interest in this field relies, among the others, on quite a remarkable mechanism: 
the amplification of quantum (incidentally vacuum) fluctuations in the metric 
of a gravitational background (see e.g. |1I6I12| '). In the presence of a change of 
regime under the cosmological evolution, the occupation number of the initial 
quantum state would get indeed amplified. As long as the change can be con- 
sidered as adiabatic, the amplification factor approaches one. Nevertheless, in 
case the change is sudden the amplification mechanism cannot be neglected. In 
such a case, even the vacuum state transforms into a multiparticle state in the 
Fock space appropriate to the new regime. The scenario clearly sounds highly 
attractive under the stronger convintion that, thanks to the great progress we 
are witnessing in the experimental application of new technologies, the ampli- 
fication mechanism may provide the possibility to detect quantum effects {e.g. 
relic gravitons) at scales considerably above the Planck one. Because of the 
semiclassical approximation underlying these studies, the natural formal arena 
turns out to be that of coherent states. These states are generated by the dis- 
placement operator D (a) (see Section 2). Squeezed states enter in the matter 
whenever the quadratic operators and ^ are involved. A squeezed state 
is generated by the action on a coherent state of the so-called squeeze operator 
defined in Section 2. As a matter of fact, in a Friedmann- Robertson- Walker 
(FRW) spacetime the behavior of matter scalar fields as well as of gravitational 
waves is governed by an equation of the time-dependent type (or time-dependent 
oscillator, TDO). Thus the problem of both particle creation and metric field 
fluctuation amplification during the cosmological evolution is reduced to that 
of solving the quantum TDO problem. The latter problem has been the sub- 
ject of several studies, mainly in connection with quantum optical arguments. 
In this paper, the main idea is to make use of the machinery built-up in 
and further developed in |I4| . The formalism would enable us to find at every 
stage information concerning the spectrum of created modes. In the case of a 
de Sitter gravitational field, the prescription straightforwardly results into full 
determination of quantities of physical relevance, such as Bogolubov coefficients 
and the phase of gravitational waves for any value of the conformal time r/. All 
these quantities are determined exactly. This has been possible by solving a non- 
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linear ordinary differential equation for an auxiliary field (j{r]) which has been 
expressed in terms of two independent solutions of the r;— dependent part of the 
D'Alembert equation for the gravitational perturbation tensor field. Another 
interesting result achieved in this paper is an estimate in terms of the auxiliary 
field a{ri) and the scale factor a{r]) of the destruction of quantum coherence due 
to gravitational time-evolution. In our approach, the use of the (nonlinear) airx- 
iliary equation in the linear parametrically excited oscillator equation for y{r]) 
(see Eq. (I2.12|) ') reveals therefore to be quite profitable and makes more strict 
and compact the formal connection between quantum optics and the theory of 
gravitational waves. The paper is organized as follows. In Section 2, after a 
description of some basic properties pertinent to the squeeze operator (|2.3|l . we 
introduce the nonlinear equation H2.9I) for the auxiliary field a{t) in terms of 
which the position and momentum operators Q and P turn out to be expressed. 
The role of the matrix element between squeezed states |a, z) of the operator 
D{a) S{z) H{t) (z) D^a) is investigated. This matrix element results to be 
evaluated in terms of the auxiliary field cr(t). It is worth noticing that in the 
expression for the matrix element three energy terms appear, one of them, for- 
mula (|2.25() . can be interpreted as the energy related to squeezed states which 
do not preserve the minimum uncertainty. In the theory of gravitational waves, 
Eq. (|2.25|) plays the role of decoherence energy of the waves. The Bogolubov 
transformation is reported whose coefficients are explicitly written in terms of 
(7. In Section 3 we discuss classical and quantum aspects of the generation of 
gravitational waves. An exact formula for the phase of a gravitational wave is 
obtained. A natural approach to the theory of gravitational waves based on 
the Kanai-Caldirola oscillator is outlined by means of an operator Hamiltonian 
expressed in terms of the auxiliary field cr(?7). In Section 4 some applications are 
displayed. Precisely, we evaluate the decoherence energy in the quasi-de Sitter 
inflationary model and standard de Sitter spacetime and the role of Bogolubov 
coefficients in terms of the auxiliary field in the particle creation mechanism is 
analyzed. Finally, in Section 5 some future perspectives are discussed. 

2 Preliminaries on the squeezed states in gen- 
eralized oscillators 

Let us recall that a squeezed state of a quantum system is defined by 

\a,z) = D{a)S{z) |0) (2.1) 

where 

D(a) = e""o-"*°o (2.2) 
is the (unitary) displacement (Weyl) operator, 

S'(z) = e5^^«'~5-*^° (2.3) 

is the (unitary) squeeze operator |2I15| . oq and a J stand for a{t)\t=to and 
a^(i)|t=to, respectively, where a{t) and a^(i) denote the annihilation and cre- 
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ation operators of the system. The complex functions a(t) and z(t) are arbitrary, 
namely 

a=\a\e"^, z = r e"t' , (2.4) 

with ip, (p arbitrary (real) c- numbers. Notice that for z = 0, Eq. (|2.1f) reproduces 
the coherent state |a,0) — D{a) |0). 
The following relations 

b ^ S^clqS ^ ao cosh r + die'''' sinhr (2.5) 
^ dlS ^ alcoahr + aoe-''!' sinhr (2.6) 

hold. This can be readily seen by applying the Baker-Campbell-Hausdorff for- 
mula. In other words, the squeeze operator S{z) induces a canonical transfor- 
mation of the annihilation and creation operators, in the sense that [6, 6^] = 1. 

2.1 Transformation of the position and momentum vari- 
ables under the squeeze operator S{z) 

From the quantum theory of generalized oscillators, it follows that position and 
momentum operators Q and P can be expressed by |13I14| 

Q=^a{a + d^), P = Vfml{^d + ea^) , (2.7) 
with d = d(t), — (t) and 

«^5.(.-f.), (2.8) 

where M — M{t) = ^, the dot means time derivative, and (the mass) m = m(t) 
is a given functions of time. The function (c-number) a(t) satisfies the nonlinear 
ordinary differential equation |18ll9j 

a + n^a^-^ , (2.9) 

(r? is specified below, see Eq. (|2.13|) '). called the auxiliary equation associated 
with the classical equation of motion 

q + Mq + uj^-{t)q^O , (2.10) 

and Lu{t) is the time dependent frequency. Via the transformation 

-i J* M{t')dt' ^^\ 

g _^ g 2 Jto y I y ^ (2.11) 

Eq. H2.10|l can be cast into the equation 

y + n\t)y^O , (2.12) 
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where 

n^{t) ^ i (alu^ - 2M - M^^ . (2.13) 

The quantum theory of the generahzed osciUator (|2.1U|) can be described by the 
Hamihonian operator jl5ll3| 

Hit) = ^ + ^mLu^Q' (2.14) 

where the canonical variables Q, P are given by Eq. (|2.7|l . 

Taking account of Eqs. H2.5|l and H2.6|l . and choosing (j) — we obtain^ 

S^QS = e''Q, (2.15) 

S^PS = 2m(j (^«^- ycr^ sinhrg + e-''P . (2.16) 

The physical meaning of the expression a — ^-a will be clarified later. At the 
present we observe only that whenever the condition 

d-^a = (2.17) 
is fulfilled, then the uncertainty product (AQ) (AP) of the variances 

i^Q) ^ ^ {Q^) - {Q)\ (AP) = ^ {P^) - {P)\ (2.18) 

attains its minimum, where the expectation value (...) is referred to coherent 
states [HI- 

We observe that the operators S'^QS and S'^PS obey the same commutation 
relation as Q and P, that is 

[S^QS, S'fPS] = [Q, P]=ih. (2.19) 

However, in contrast to what happens for the operators b and 6^ (see Eqs. 
(E^J-lESJl), the operators S'^QS and S'lPS are not Hermitian conjugate. We 
point out that the property of Hermitian conjugation is enjoyed by the operators 
S^QS and S^PS in the case in which the condition H2.17|l is valid. 

A possible physical interpretation of the properties H2.15|l and (|2.16|l is the 
following. For a quantum system governed by a Hamiltonian preserving the 
minimum wave packet (i.e., the condition (|2.17|l holds) Eqs. H2.15|l and H2.16|l 
become S^QS — e^Q and S^PS — e~^P, respectively. If 1-0) is the state of the 
system under consideration, then — S{r) lip) represents the same system 

^Stoler 1^ saw that, generally, states of the type \a, z) do not describe wave packets relative 
to the minimum value of the product (AQ) (AP), where A means the variance operation (see 
Eq. trr^ ). The state \a, z) can describe a wave packet of minimum uncertainty only if z 
is real (0 = 0). In the framework of quantum generalized oscillators, this corresponds to the 
condition I2.17i . 
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squeezed in the space of the position Q by a factor e and expanded in the 
space of the momentum P by the factor e''. In fact, we deduce 

e-'' (V'' I Q I = (V' I Q I ^> , e^{^'\P\^') ^{^\P\i^) . (2.20) 

Now, we shall evaluate a matrix element involving the operator H{t) (see 
Eq. H2.14|l ) in the context of squeezing of a quantum system. In doing so, let 
us consider the following expectation value between squeezed states: 

{a, z I D(a) S{z) H(t) S\z) D\oi) | a, z) = 

= (0 I S\z) D\a)D{a) S{z) H{t) S\z) D\a)D{a) S{z) \ O) = 



= (0 I H{t) I 0) - ^ (0 I P2 I 0) + \mu;\t) (O | | O) 

(2.21) 

where Eqs. (|2.1I) and (|2.14|l have been employed. 

The expectation values on the right hand side of Eq. I|2.21|l can be evaluated 
from Eq. ifTTji . They read 



(0|Q2|o)^la2, 



hm 



1 

4^ 



M 



(2.22) 



(2.23) 



where ^ and a are described by Eqs. (|2.8|) and H2.9|l . With the help of Eqs. 
and Eq. itT^ takes the form 

(a, z I D{a) S{z) H{t) S\z) D\a) | a, z) = (0 | H{t) \ 0) = 



h h 
8^^ 2 



M 



(2.24) 



The term in the square bracket corresponds to the vacuum expectation value of 
the kinetic energy of the system, while the last term is related to the vacuum 
expectation value of the potential energy. 
The quantity 



Enm 4- 



h ( M 



(2.25) 



(NM=non- minimum) can be interpreted as the energy associated with the squeezed 
states which do not satisfy the criterium of minimum uncertainty [Enm 0). 
When the criterium is verified, then E^m = 0. In such a case Eq. H2.24|l can 
be written as 



(a, z I D{q) S{z) H{t) S''{z) D\a) | a, z) = (0 | H{t) \ 0) 



h 
8^ 



(2.26) 
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Hence in the minimum uncertainty situation the vacuum expectation values of 
the kinetic energy and the potential energy turn out to be proportional to 
and respectively. It is noteworthy that, in general, all the energies appearing 
in Eq. (|2.21|l can be expressed in terms of the auxiliary field ait) obeying the 
auxiliary equation H2.9|l . This enhances the convenience of our approach to the 
study of cosmological quantum effects based on the theory of Eq. (|2.9|) , which 
we are going to develop in Section 3. 

We remark that the quantity ("y-o" — cr) is connected with the expectation 
value of the operator {Q, P} = QP + PQ between vacuum states, i.e. 

(0|QP + PQ|0) =2no-(^ycr-(T^ . (2.27) 

Thus the minimum uncertainty requirement ^-a = a implies that the expecta- 
tion value (0 \ QP + PQ \ 0) is vanishing. 

2.2 A link between Eq. (j2.12|) and the auxiliary equation 

For later convenience (see Section 3), we shall report a result establishing a 
relationship involving the solutions of the (linear) equation of motion and the 
(nonlinear) auxiliary equation 



V?{t)G^^ , (2.28) 



where k is a constant. If y\ and 7/2 are two independent solutions of Eq. (|2.12|) . 
then the general solution of the auxiliary equation H2.28|l can be written as 

a=(Ay2 + By2 + 2Cyiy2)5 , (2.29) 
with A, B, C arbitrary constants such that 



^B-C^^^ (2.30) 

where Wq = Wo(yi, 2/2) = 2/12/2 ~ 2/12/2 = const is the Wronskian. 

It is worth remarking that from the theory of the auxiliary equation H2.28|l 
a phase can be given by the real function 0{t), 

m^lJ^ . (2.31, 

(See |18I19I20| : for some applications: |16ll3j .') Here we shall suggest the pro- 
cedure which can be used to compute the above integral in general cases. To 
this aim it is convenient to introduce the function 

m = ^e''"yi{t)-^e'^' y2{t) , (2.32) 

where a,/3 are real numbers and 2/1,2/2 are the two independent solutions ap- 
pearing in Eq. (|2.29|l . We have 

\^P{t)\^ ^Ayl+Byl- 2^AB cos^o 2/i2/2 , (2.33) 
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where 6*0 = a - p. Comparing Eq. (|2.33|) with Eq. H2.29|l we get C 
— \JAB cos 610, so that the condition H2.30|l becomes 



ABsin2^?o = -^, (2.34) 



from which k 7^ whenever sin^o 7^ 0. Hence, the auxihary field a can be 
expressed by 

a\t) = 1^(01' . (2.35) 

To calculate the phase d(t) corresponding to the solution (|2.29() of Eq. (|2.28() . 
we look for a function F(t) defined by 

so that 

The numerator in Eq. H2.37() can be elaborated to give 

V^V-* - V'V'* = 2i\/ABsin6ioTyo 2iVK (2.38) 

where Eq. has been exploited. Substitution from Eq. in Eq. l|07jl 
thus yields 

F = . (2.39) 
Then the phase B(t) is determined by integrating (|2.31|) . namely 

/■* dt' 

= ^ = m - F{t^) , (2.40) 
where F(t) is provided by Eq. H2.36|) . We shall recall this general result later. 

2.3 The Bogolubov coefRcients in terms of a 

By resorting to the operators 

where oq = a(io) is the (mode) annihilation operator in the Schrodinger rep- 
resentation, combining Eqs. 12.41|l and H2.7|l we can derive the Bogolubov 
transformation 

a{t) = pi{t) ao + v{t) 4 (2.42) 
whose coefficients are expressed by 



f^W = \h [-< + crl , v{t) = J- [-1^ crl , 

V ZlOqWIq V TO / V ^WqTOo V TO / 

(2.43) 
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with mo — m{to), ojq — o-!{to), and ^ given by Eq. 1)2. 8|l . We note that the Bo- 
golubov transformation can be naturaUy embedded into the relations (|2.5() and 
H2.6|l . where the operators b, can be identified with a{t) and a^{t). Equation 
entails 

Im|'-|H' = 1- (2.44) 
On the other hand, the uncertainty product can be formulated as follows |14| : 

(AQ) (AP) = ^^l + 4a2 = 

= ^Kt)-Ki)l |M^) + K^)|> J. (2.45) 

The uncertainty formula H2.45|l is closely related to the concept of coherent 
states for the generalized oscillators. Such coherent states were constructed by 
Hartley and Ray in 1982 [2] taking account of the Lewis- Riesenfield theory |16| . 
These states share all the features of the coherent states of the conventional 
(time-independent) oscillator except that of the uncertainty formula, in the 
sense that the product (AQ) (AP) turns out to be not minimum. A few years 
later, Pedrosa showed that the coherent states devised by Hartley and Ray are 
equivalent to squeezed states 

3 Classical view and quantum theory generation 
of gravitational waves via the Kanai-Caldirola 
oscillator 

In this section we shall develop a model of propagation of gravitational waves 
based on the application of the auxiliary equation H2.9() for the function cr(ry), 
in terms of which the Bogolubov coefficients can be built-up. The Bogolubov 
transformation is a basic concept in the theory of particle creation in external 
fields. The created particles do exist in squeezed quantum states According 
to relic gravitons created form zero-point quantum fluctuations during cos- 
mological evolution should now be in strongly squeezed states. In this context 
the generation of gravitational waves is of fundamental importance. 

The theory of generation of gravitational waves in the inflationary universe 
scenario is based on the action [11] 

^-i6b/^(^)^^'^' (^-1) 

where f{R) is an arbitrary function of the scalar curvature R. The theory 
defined by the above action is conformally equivalent to a pure Einstein theory 
with scalar-field matter. In linear theory, the gravitational waves decouple from 
the matter field, so that the main problem is to fix the background model and 
to desume the relation between the conformal metric. 
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Starting from Eq. H3.1|l . by varying the action with respect to the gravita- 
tional perturbation field hj, the equation of motion 



h" + 2—h' + {2K-A)h^0 (3.2) 
a 

is obtained, where prime denotes a(rf) — \J'§^o.{r]), h = h{ri,^) is each 

component of /i* , A stands for the Laplace-Beltrami operator, and K means 
the space curvature. By separating in h{ri, ~x) the dependence on 77 from the 
dependence on ~x, we can write h ^ ho{r]) hi{~x), so that 

[A + {n^~K)] hi{lc)^0, (3.3) 

and ^ 

h'^ + 2^h'Q + {n^ +K)ho^O. (3.4) 

Equation (|3.4|l can be applied to describe the evolution of gravitational waves 
in any state of the evolution of the universe, even when resorting to higher 
derivative theories of gravity . The elimination of the first derivative in Eq. 
H3.4[l leads to the equation (called master equation in ^) 

y" + [{n^ + K)~Vifj)]y^O , (3.5) 

with 

a a{r]o) 

We remark that Eq. (|3.4|l can be regarded as the equation of motion of an 
oscillator with time-dependent mass m and constant frequency cj, 

q" + !!^q' + Cj\^0 , (3.6) 
m 

which is described by the Hamiltonian 

The quantum theory of gravitational waves is therefore equivalent to the quan- 
tum theory of the Kanai-Caldirola oscillator |23l24j . The formal analogy is 
realized upon the identification: 

,2 -2 _ 2 



a , Lu = n + K , q '--^ ho ■ (3.i 



3.1 Exact solution of the parametrically excited oscillator 
and its associated auxiliary equation 

We are mainly interested in the period under which the universe accelerates, 
namely in its inflationary stage. Infiation is defined to be a period of accelerating 
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expansion. During such a stage, the universe expands adiabatically and the 
Friedmann equations can be exploited ^Jj. The prototype of the models of 
inflationary cosmology is based on the de Sitter spacetime, which is a very 
interesting case concerned with the constant Hubble rate and the scale factor 
given by 

<V) = —^ , (3.9) 
where 77 < and Hq denotes the Hubble constant. In the de Sitter case, Eq. 



takes the form of Eq. (I^H)) . with n'^{'q) = 
solution 



2 

7p 



It admits the general 

y = \/— n?7 KiJ3{—nrj) + K2J_i{—nf]) , (3.10) 

where J3 , J_ 3 are Bessel functions of the first kind and the arbitrary constants 
Ki and K2 are determined once the initial conditions are imposed. Then from 
Eq. H2.29(l we infer that Eq. (|2.28() , which now reads 



(3.11) 



(3.12) 



is exactly solved by 

(T = \/"^my AJi{—nrj) + BJ^3{—m]) + 2CJi(— n?7) J_i{—nrf) 
where constants A,B,C satisfy the condition l|2.30|l . 



3.1.1 The phase of the auxiliary field <j{t) 

As we have shown in Subsection 2.2, to calculate the phase 9 corresponding to 
the solution (|3.12(l of Eq. (|3.11(l it is convenient to introduce the function 



(3.13) 



where a,/3 are real numbers. Hence the auxiliary field can be expressed via 
(T^ = I -01^ provided that 

ABsin^Oo^-^, (3.14) 



where 6*0 = a — (3 and Wq = Wo{y/—nrjJ^ {—nrj), ^/—nr/J^ {—mf) ). Therefore 
it is an easy matter to see that in the case under consideration the phase 9 
determined by integrating Eq. (|2.4UI) is given by (z = — n?]) 



0(77) = 



dz' 



where the function F{z) is defined by 



-^-n[F{z)-F{z,)] 



F{z) = 



1 



:ln 



(3.15) 



(3.16) 



and ip is given by Eq. H3.13|l . 
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3.2 On the phase of gravitational waves 



A deep discussion on the phase of gravitational waves is contained in j25| , where 
this topic is dwelt on both at the classical and quantum level. 

Here we confine ourselves to tackle the problem classically. The study of the 
phase of gravitational waves by a quantum point of view will be done elsewhere. 

For our purpose, first we observe that the Bessel functions J 3 and J_3 can 
be explicitly written as follows j^Hl 



J3(Z) 

Then, the auxiliary field 



2 / . cosz^ 
sm z H 

TTZ \ z 



A cos z 



can be written as 
2 



/ cos z \ 2 

+B[sinz+ ^—j ~2C 



sm z 



( . COSZN 
[smz+^j 



(3.17) 



where t/i = ^/zJa, 1/2 — \fzJ_3 and the condition H2.30|l is understood. 

Now we shall see that the phase Qgw of the primordial gravitational waves 
can be obtained by Eq. H3.17|l under the choice 



A = S7^0, C = 



(3.18) 



and assuming asymptotically negative values of the conformal time. In doing 
so, Eqs. IpTT^ . (pTTTjl provide 



2A 

TT 



2A 



with 



A^ 



2n 



(3.19) 



(3.20) 



In the case -n? ^ l^l??)!: the high frequency waves, e.g. the solutions of the 
equation 

y" + n^y = , (3.21) 
correspond to the following behavior of the gravitational perturbation field h: 



h{rj) — — sin {nq + p) 
a 



(3.22) 



(p is an arbitrary phase). In an expanding universe, the amplitude h of the waves 
decreases adiabatically for all 77. The result represented by the calculation of 
the phase (|3.19l) of relic gravitational waves suggests one to interpret formula 
H3.15|l as the phase of gravitational waves not only in the case of the primordial 
cosmological scenario formally corresponding to 77 ^ —00. Anyway, this subject, 
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which goes beyond the scope of the present paper, deserves further investigation. 
Here we recall only that the general solution of Eq. (|2.12() , which holds in the 
case of the de Sitter cosmological model, can also be written as 

V = -5=fT(?y) cos [V^e{ri) + 6] (3.23) 
v K 

where 9 is given by Eq. 13.15|l . c is a Noether invariant of Eq. 12.121) . and 5 is 
an arbitrary constant. Hence for any conformal time "q in the interval (— oo, 0), 
the amplitude h of the gravitational perturbation field can be expressed by 

h ~ —criv) sin [y/K6{r]) + const] . 

3.3 Quantum gravitational waves: Theory in terms of the 
auxiliary field 17(77) 

The approach to the study of gravitational waves we present in this paper is 
developed starting from the classical Kanai-Caldirola Hamiltonian j23l24| . Eq. 
H3.7|l . As we have already pointed out, Eq. (|3.4|l can be regarded, in fact, as 
the equation of motion of an oscillator with time-dependent mass m = and 
constant frequency u) = + K (see Eq. (I3.()|l ) described by the Hamiltonian 
1)3. 7|l . So, we make the fundamental identification of the whole temporal part ho 
of the metric fluctuation amplitude as the basic " coordinate" variable to quan- 
tize as such. (Recall that our procedure for the quantization of gravitational 
waves is based on the identification of y — aho as the variable to quantize). 
As a consequence, the quantum theory of gravitational waves turns out to be 
completely equivalent to the quantum theory of the Kanai-Caldirola oscillator 
which can be described by the quantum version of the Hamiltonian (|3.7|) . On 
the grounds of what we learned in Section 2, the above identification suggests 
a route which can be successfully pursued whenever we are interested in the 
characterization of physical effects (quantum decoherence, squeezing, particle 
production etc.) emerging from the study of inflationary models in the early 
universe. Section 4 will be devoted to a preliminary exploration of the effective- 
ness of the idea in the context of expanding universe cosmological models. 

4 Applications 

By taking full advantage of the formalism introduced in Section 2, we are in the 
position to study the dynamical system of the cosmological interest which are 
described by time-dependent oscillators. In doing so, a key point is the char- 
acterization of constants A, B, C in Eq. (|2.29|) . It is concerned with the initial 
(and boundary) conditions. All dynamical aspects under time evolution are en- 
closed into the function cr, which obeys the second order nonlinear differential 
equation (|2.9() . A general condition on A,B,C is provided by Eq. H2.30|l . It 
is not enough, however. Specification of the value and the first time derivative 
of (7 at fixed time is thereby needed. Another condition is associated with the 
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requirement that at the initial time rji the time-dependent annihilation and cre- 
ation operators, derived from Eq. 12.7II . go into the standard Dirac-like form, 
Eq. H2.41|l . As for the final condition, it is helpful to reveal that in most cases 
we want the state at initial time to correspond to a vacuum state. This can be 
achieved easily under the minimization requirement for Ej^m — f i,^^'^ ~ ^) ■ 
Indeed, since it provides a measure of the decoherence at the time 77, it has to 
be vanishing when referring to a vacuum state at the initial time rj = rji. Under 
these circumstances, the whole set of initial conditions for a is given by 

a{to)-{AuY'^'^0, (4.1) 
a{to) - ^a(io) = . 

In a cosmological framework of the FRW type, above system is translated into 

In order to proceed with concrete analysis, it is very customary to resort to 
the spatially flat inflationary model based on the de Sitter metric. However, 
a more general and realistic description of the inflation may be provided by a 
quasi-de Sitter spacetime (see e.g. |2H1)- In this case, the Hubble rate is not 
exactly constant but, rather, it weakly conformal changes with time according 
to H' — —ea^H"^ (that is, aa" = (2 — e) a'^H'^ = (2 — e) a'^) where e is a constant 
parameter. When e vanishes one gets just the ordinary de Sitter spacetime. For 
small values of e, a quasi-de Sitter spacetime is associated with the scale factor 

^^n) = . (4.3) 

H{1 - e)r) 

(ry < 0). In the quasi-de Sitter spatially flat scenario, Eq. ()2.12f) reads 

2 (2 + 3e) 



y 



n 



y = (4.4) 



and can be solved in terms of Bessel functions. Precisely, one has the two 
independent solutions 



2/1 = V-myJ^ (-nri) , y2 = V^rujY^ {-nrj) (4.5) 



where ^ = y j ~^ ^(t-i)^^ ■ procedure outlined in the previous sections can 
be applied and we are led to the introduction of the basic function 

= {AJl {-~nri) + BY^ {-n-q) + 2CJ^ {^nr]) {-nr]) } ^ (4.6) 
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where A,B,C are determined by means of the system (14.211 . r]i denoting the 
conformal time of the beginning of the inflation. Once we are interested in a 
situation in which the system started very far in the past in a vacuum state, 
the Bessel function expansions 



Jv{~nri) 



2 

nnrj 



cosi— n?7 — — TT -) + O 

2 4 



1 

nr] 



Yui-nrj) - \ 

" nnr] 



sm 



-nrj- -TT- 



TT, / 1 

T +0 — 

4' \nr] 



for v fixed and nr] — > —oo assists us in finding suitable constants A, B, C. 
By taking arbitrary asymptotically negative initial times, the leading terms 
of Bessel functions J^, give rise to the following behavior for the function: 



a{r]) 



- {A + (B - A) sin^f-njy - -tt - -) + Csin(-2?7?7 - i/tt - -) + O ( — 
_ I o /f ^ nr] 



Once the limit r]i << is concerned a natural choice is given hy A = B = 
C = (recall that we already found this result for the case v — 3/2 associated 
with the standard de Sitter metric background). So we obtain 



'^(^) = V"!'? {4{-nr]) + Y;^{-nr])y 



-^V \Hl{~nr])\ 



(4.7) 



In the light of our previous results, the decoherence energy Ej^m at the time 77 
of gravitational waves in a quasi-de Sitter model of inflation can be evaluated 
by inserting (|4.7|) into formula H2.25|l . It then results 



Enm — 



a{r]) 



1 2 



(l-e)77. 



where 



{l-e)r] 
Moreover, since 




2 \Hl{-nr])\ 



3-e 
2(1^ 



\Hl{-nr])\ 



(4.8) 






r , «' 1 




— i 


a cr 


} 




a 






\ ' °' 1 




— i 


a cr 


} 




a 





(4.9) 
(4.10) 
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at an arbitrary time 77 the Bogolubov coefficients are given by 



2n 




1 

2^ 



1 

2^ 



(l-e)77 



(4.11) 
(4.12) 



with a and <t' + ^ furnished by Eqs. 14.7|l and 14.8|l . respectively. Finally, 
the phase 9 can be evaluated. Due to Eq. (|2.4Q(I . we get 



i(ry) 



'i In 



■tjj{r]) 



(4.13) 



where 



That is, 



3I I'' 

l»7i 



where 9]^ denotes the phase of the Hankel function Hi . 

It is now instructive to focus on a standard de Sitter inflation. In this case 
e = 0, = I and 

1 



/2n 



1 



1 



So, the integration of ^ is straight, i.e. 



= An' 



2 2 



1 



drj — An 



tanh ^ (nrj) 



(4.14) 



(4.15) 



up to a constant of integration. The above two formulas provide the standard 
exact (and normalized) solution for the quantum fluctuations of a generic mass- 
less scalar field during a de Sitter inflation. Yet, it is interesting to consider 
a very long inflation by letting the system evolve towards 77 —> 0. In this case 
Enm simply reads 



Enm{vi) 



M 



4 1 + 7^^77^ 



(4.16) 



Interestingly, the decoherence energy at 77 = is finite. Moreover, in the stan- 
dard de Sitter phase, from Eq. H4.14|l we obtain the Bogolubov coefficients 



1 \J n'^Vi 

2 y^l + n^jf' 



2 + 772772 



1 + 


r;2 1 


i 

nrj 


1 - 




i 

nrj 



(4.17) 
(4.18) 
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{p?rif ^ 1) which in turns imply that 

4 V ^7 V^J 

particles are created out the vacuum at the time rj. 

In more refined studies of cosmological effects in the expanding Universe, 
it turns out to be useful to introduce a cosmological model which allows one 
to take into account different evolutionary phases of the universe. Once the 
model has been specified and the equation H2.12(l solved, one can get an insight 
into physical effects associated with different cosmological stages. For instance, 
one can consider a simple cosmological model which includes the inflationary 
(i), radiation-dominated (e) and matter-dominated (m) epochs |29| . The scale 
factor has the following dependence on the conformal time: 

< aeirj) = T^^^ , 1le<1l<Vm , (4.20) 

where Hq denotes the Hubble constant at the inflationary stage, and rji repre- 
sents the beginning of the expansion. In order to determine the ?7-dependent 
amplitude /i„ for each epoch we have to solve Eq. H2.12|l with the corresponding 
varying frequencies, namely 

/ 2 2 V 

Q,i ^ In - — j , fie = « , ^rn = 

Notice that the frequency Qe is constant, while fli and flm are varying with 
the same temporal dependence. The frequency fli characterizes the de Sitter 
era. The related equation of motion H2.12(l has already been solved. Taking care 
about matching data at rje and ry™ is needed for the knowledge of the complete 
form of a. The associated cr's and their derivatives have to join continuosly at 
?7e and 77™, in fact. This step is needed to obtain all the physical information 
implied in formulas for the Bogolubov coefficients, the decoherence energy, the 
gravitational phase, and squeezing. Having in mind our previous discussion, 
employing the model by considering a quasi-de Sitter phase is straight. 

In general, the vacuum expectation value of the number operator and the 
other quantities of the physical interest vary slowly with time if the expansion 
rate becomes arbitrarily slow. In case the espansion is stopped one should be 
able to recover time-independent Dirac operators. However, the circumstance 
does not mean that Bogolubov coefficients trivialize. This is because loss of co- 
herence previously occurred due to the expansion dynamics. A typical situation 
may be that a goes to zero but a does not. Enm goes to zero, indicating that 
when expansion is stopped the time-dependent gravitational pumping stops as 
well and there is no further decoherence. If expansion stops from time rja to 



iv + rim - ^rieY 



(4.21) 
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time 77h, then V77 G [rja, rji,] one gets 




2na?{rii) \2a{r]a) 



which impUes 



n— — Ta(r\a) 



1 



(4.22) 
(4.23) 

(4.24) 



To get a more clear insight into the resuhs achieved in this section, a few 
comments are in order. Specifically, our formula for the decoherence energy 
associated with the dynamical evolution of the gravitational fluctuation modes 
on a background of the FRW type, is expressed in a very compact form in terms 
of the auxiliary field a-{r]) and the scale factor a{ri). On the other hand, the de- 
coherence energy plays an essential role in the relationships for the Bogolubov 
coefficients. This aspect makes explicit how the energy lost owing to the deco- 
herence effect may be exploited to excite the vacuum state of the model under 
consideration. Moreover, in our framework this mechanism would be quantified 
in a compact way by means of the formula 



1 

2^ 



2 
n 



(4.25) 



where E^m is the decoherence energy (to — a?). With respect to other works, in 
our paper the role of the decoherence energy is made manifest. Furthermore, we 
observe that, remarkably, the auxiliary field a is nothing but the time-dependent 
amplitude of the mode solutions to Eq. (|2.12|) for the redshifted gravitational 
field fiuctuations. 



5 Concluding remarks 

The main results achieved in this paper have been presented and widely dis- 
cussed in the Introduction. Therefore we shall conclude by making some final 
comments concerning challenging perspectives which should be dwelt upon in fu- 
ture developments. The evolution equation of a mode with comoving wavenum- 
ber n reduces to the harmonic oscillator equation with time-dependent mass and 
constant frequency. The approximation behind computations leading to the re- 
sult actually are applicable only to the infrared region. On general grounds, one 
therefore expects that predictions for observables may depend sensitively on the 
physics on the length scales smaller than the Planck one. In order to take into 
account trans-Planckian physics, it has been recently suggested to make use of 
effective dispersion relations (see e.g. [2Z|)- The linear dispersion relation is 
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thus replaced by a nonlinear one, n^^y = a? {rf) F"^ [n / a) , where F{n/a) is an 
arbitrary function required to behave linearly whenever n/a (= k) is below a 
certain threshold. A time dependent dispersion relation thus enters in the mat- 
ter. As a consequence, the underlying dynamical model turns out to be that 
of the harmonic oscillator with both the mass m {=a?) and the frequency lo 
(= TLeff) depending on time. In principle, its quantization can still be pursued 
by resorting to the formalism of Section 2 and it will be studied in detail else- 
where. Nevertheless, a comment is in order. In the light of the discussion in 
|14) . one might wonder, in fact, on whether or not in the cosmological frame- 
work the minimum uncertainty criterium can be satisfied under time evolution 
for some physically reasonable function F. It is straightforwardly seen that this 
is not the case, generally speaking. Indeed, in the cosmological framework the 
criterium reads as a^Ueff = const and implies a purely cubic function F, say 
F = a^n^ /a? . As a consequence, the uncertainty relation can be minimized only 
approximately. It is worth noting that this happens in the large wavenumbers 
limit of a special case of the generalized Corley-Jacobson dispersion relation 
introduced in ^ (see Eq. (22) in |271). 
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